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ABSTRACT 


The  successful  application  of  holography  to  the  study  of  three  dimensional 
flow  fields  due  to  phase  objects  has  been  reported  in  the  literature.   The 
present  report  extends  this  technique  to  the  study  of  density  fields  around 
opaque  bodies  as  would  normally  be  encountered  in  wind  tunnel  experiments. 
The  density  field  around  a  10  degree  half- angle  cone  at  zero  and  ten  degree 
angle  of  attack  has  been  investigated  by  means  of  the  finite  fringe  holographic 
interferometry .   The  three  dimensional  density  field  obtained  from  the  reduc- 
tion of  the  interferograms  was  found  to  agree  with  that  obtained  from  an 
analytical  solution  of  the  governing  equations. 
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INTRODUCTION 

This  is  a  report  on  holographic  investigations  directed  towards  the 
determination  of  complex  three-dimensional  density  fields.   The  report 
consists  of  essentially  Appendix  I  which  includes  the  thesis  of  R.  C. 
Jagota  on  the  determination  of  density  fields  about  a  cone  at  angle  of 
attack  in  a  wind  tunnel. 

A  paper  based  on  this  thesis,  "Finite  Fringe  Holographic  Inter- 
ferometry  Applied  to  a  Right  Circular  Cone  at  Angle  of  Attack,  by  R.  C 
Jagota  and  D.  J.  Collins,  has  been  accepted  for  publication  by  the  Journal 
of  Applied  Mechanics. 

Further  work  on  the  determination  of  the  density  in  corner  flows  was 
begun  towards  the  end  of  the  contract  year.   This  work  is  directed  at  determining 
the  flow  in  a  wing  fuselage  junction  at  transonic  speeds. 

A  preliminary  program  using  a  flat  plate- fin  in  the  Naval  Postgraduate 
School  wind  tunnel  has  been  started  at  Monterey.   This  preliminary  program 
is  designed  to  solve  some  of  the  holographic  problems  associated  with  the 
transonic  test  that  are  to  be  conducted  at  NSRDC.   Further  reports  on  this 
work  will  be  forthcoming  in  December. 
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ABSTRACT 

The  successful  application  of  holography  to  the  study  of 
three  dimensional  flow  fields  due  to  phase  objects  has  been 
reported  in  the  literature.   The  present  report  extends  this 
technique  to  the  study  of  density  fields  around  opaque  bodies 
as  would  normally  be  encountered  in  wind  tunnel  experiments. 
The  density  field  around  a  10  degree  half-angle  cone  at  zero 
and  ten  degree  angle  of  attack  has  been  investigated  in  the 
Naval  Postgraduate  School  supersonic  wind  tunnel.   In  these 
experiments,  the  finite  fringe  method  for  the  production  of 
interferograms  has,  for  the  first  time,  been  applied  to 
holographic  interf erometry .   The  density  field  obtained  from 
the  reduction  of  the  interferograms  was  found  to  agree  with 
that  obtained  from  an  analytical  solution  of  the  governing 
equations  as  reported  by  D.  J.  Jones  in  Reference  1.   The 
computer  program  used  for  the  reduction  of  the  interferograms 
has  been  evaluated  for  the  effect  of  the  presence  of  the  cone 
and  the  shock  wave  and  has  been  found  to  yield  stable  and 
accurate  results. 
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I.   INTRODUCTION 

For  over  three  decades,  the  Mach-Zehnder  interferometer 
has  been  widely  used  and  has  come  to  be  accepted  as  a  stan- 
dard tool  for  the  measurement  of  the  flow  field  around 
bodies.   Until  the  advent  of  the  laser,  its  use  required 
extreme  stability  from  vibrations  and  precision  optical 
components.   Even  with  a  laser  as  a  light  source,  the 
resulting  interferograms  are  purely  two-dimensional  in 
nature  in  that  they  provide  fringe  information  for  only  one 
direction  of  view  through  the  wind  tunnel.   Furthermore, 
due  to  the  limit  imposed  by  existing  camera  shutter  speeds, 
the  interferogram  obtained  does  not  truly  represent  the 
instantaneous  flow  field  around  the  body  but  rather  is  an 
average  over  the  interval  of  shutter  opening. 

The  advent  of  holography  together  with  the  development 
of  the  high  power  Q-switched  laser  has  introduced  a  new 
dimension  in  the  field  of  flow  measurement.   Apart  from 
removing  the  need  for  precision  optical  quality  components 
in  the  case  of  double  exposed  holograms,  the  use  of  holog- 
raphy has,  for  the  first  time,  enabled  a  three-dimensional 
view  of  the  flow  field  to  be  obtained  from  a  single  hologram 
The  high  power  of  the  laser  and  the  Q-switch  mechanism  have 
together  enabled  film  exposure  times  in  the  region  of  20 
nano-seconds  to  be  attained  thus,  in  effect,  "freezing"  the 
flow  during  the  hologram  production  process . 


Techniques  for  the  application  of  holography  to  inter- 
ferometry  have  been  reported  in  Reference  1  (Heflinger 
et  al,  1966)  and  in  Reference  2  (Brooks  et  al,  1966). 
Holographic  interferometry  has  been  used  to  determine  the 
three-dimensional  asymmetric  field  produced  by  free-jet  flow 
(Reference  3).   In  this  report,  a  series  of  holograms  placed 
in  a  90  degree  sector  around  the  jet  were  obtained  simul- 
taneously from  a  single  exposure.   Interferograms  were  then 
obtained  from  the  holograms  at  suitable  angular  orientations 
within  a  90  degree  field  of  view.   This  provided  sufficient 
information  to  determine  the  entire  flow  field  since  the  jet 
was  planar  symmetric.   Reduction  of  the  fringe  data  to  den- 
sity information  was  accomplished  by  expressing  the  density 
in  a  series  of  orthogonal  polynomials  and  Inverting  the 
fringe-density  integral  equation  by  a  computer  program  using 
a  scheme  reported  in  References  4  and  5. 

In  the  present  report,  the  applicability  of  holography 
interferometry  has  been  extended  to  those  three-dimensional 
density  fields  produced  by  opaque  bodies  as  would  be 
encountered  in  wind  tunnel  experiments.   The  experiments 
were  conducted  in  the  Naval  Postgraduate  School  supersonic 
blowdown  wind  tunnel  at  a  Mach  number  of  2.8  using  a  ten 
degree  half-angle  cone  at  zero  and  ten  degrees  angle  of 
attack.   In  this  case,  the  multiple  hologram  technique  was 
not  employed.   Since  the  flow  fields  investigated  were  steady 
the  procedure  adopted  was  to  obtain  holograms  to  cover  var- 
ious sectors  of  the  flow  field  by  rotating  the  cone  within 


the  wind  tunnel,  this  being  a  method  that  would  be  fully 
applicable  to  any  conventional  wind  tunnel  setup. 

Whereas  in  an  infinite  fringe  field,  the  fringe  infor- 
mation is  obtained  in  a  discontinuous  manner,  that  is,  at 
the  points  of  maxima  and  minima,  the  interferogram  from  a 
finite  fringe  field  can  be  evaluated  at  any  desired  point 
quite  readily.   This  is  of  particular  significance  in  the 
case  of  density  fields  where  the  variations  in  density  are 
small.   In  the  present  case,  the  total  variation  in  fringe 
number  was  generally  in  the  region  of  only  one  to  two 
fringes.   The  finite  fringe  method  was  also  selected  since 
it  is  less  sensitive  to  vibrations  than  is  the  infinite 
fringe  field,  a  problem  encountered  in  Reference  3-   The 
finite  fringe  field  was  produced  by  a  translation  of  the 
diffusing  screen  in  the  scene  beam  a  distance  of  between 
.025"  and  .030"  between  the  two  exposures  of  the  hologram. 
A  vertical  fringe  field,  obtained  by  horizontal  translation 
of  the  diffuser,  was  chosen  in  spite  of  the  greater  effort 
required  to  reduce  the  interferograms  at  a  particular  section 
owing  to  the  greater  uniformity  and  clarity  of  the  fringes 
thus  obtained  as  compared  to  horizontal  fringes.   Fringe 
data  obtained  from  the  interferograms  was  reduced  to  density 
using  the  same  computer  program  as  used  in  Reference  3.   The 
density  field  obtained  experimentally  was  compared  with  the 
analytical  solution  reported  in  Reference  6,  interpolating 
the  results  therein  for  the  actual  mach  number  of  the  experi- 
ments.  Good  correlation  was  achieved,  thus  demonstrating 


the  utility  of  holographic  interferometry  for  normal  wind 
tunnel  experimentation. 


II.   EXPERIMENTAL  APPARATUS 

A.  THE  WIND  TUNNEL 

This  investigation  was  conducted  in  the  Naval  Postgraduate 
School  blow-down  supersonic  wind  tunnel.   The  wind  tunnel  is 
a  blowdown-to-atmosphere  facility  which  has  a  test  section 
H    inches  by  4  inches  in  cross-section  and  6  inches  in  length. 
Interchangeable  nozzles  are  available  for  producing  nominal 
test  section  mach  numbers  of  4.0,  2.8,  and  1.7.   The  nominal 
run  time  is  5  minutes  at  a  mach  number  of  2.8  with  a  maximum 
stagnation  pressure  of  about  105  pounds  per  square  inch. 
For  this  series  of  experiments,  the  side  walls  of  the  wind 
tunnel  were  replaced  by  plexiglas  walls  two  inches  thick 
having  a  refractive  Index  of  1.49  so  that  a  complete  view  of 
the  flow  in  the  test  section  could  be  obtained. 

B.  THE  EXPERIMENTAL  SETUP  FOR  PRESSURE  MEASUREMENT 
Figure  1  is  a  schematic  representation  of  the  equipment 

arrangement  used  for  pressure  measurement.   The  static 
pressures  on  16  out  of  a  total  of  30  orifices  on  the  surface 
of  the  cone  were  measured  by  a  Scanivalve  (Model  48J4-10^5 
with  a  0-10  psia.  transducer)  which  was  calibrated  in  inches 
of  mercury  vacuum.   The  measurements  were  recorded  on  a 
Honeywell  2106  twelve  channel  Viscicorder.   The  scanivalve 
was  capable  of  scanning  15  ports  per  second.   Photographs  of 
the  experimental  setup  are  provided  in  Figures  2  and  3-   The 
plexiglas  side  walls  of  the  wind  tunnel  with  the  Scanivalve 


mounted  on  top  of  the  test  section  are  shown  in  Figure  2. 
The  power  supply,  control  box  and  Viscicorder  can  be  seen  in 
Figure  3. 

C.   THE  HOLOGRAPHIC  ARRANGEMENT 

A  schematic  of  the  holographic  arrangement  employed  is 
shown  in  Figure  4  and  photographs  thereof  are  in  Figures 
5(a)  and  5(b).   The  equipment  was  mounted  on  a  table  that 
rested  on  a  portion  of  the  floor  that  was  vibrationally 
isolated  from  the  rest  of  the  building.   The  monochromatic 
light  source  used  was  a  Korad  K-l  pulsed  ruby  laser  operating 
at  a  wavelength  of  69^3  Angstroms  together  with  a  Pockels 
cell  Q-switching  device.   The  resulting  effective  exposure 
time  was  about  20  nanoseconds  and,  with  a  cavity  length  of 
the  laser  of  37  cms.,  a  coherence  length  of  about  2  cms. 
was  obtained.   A  side  band  geometry  was  employed  to  obtain 
the  holograms  of  the  flow  field  around  the  cone.   A  diffuse 
glass  in  the  path  of  the  scene  beam  was  used  to  produce 
light  field  holograms,  which  because  of  the  diffusion  of 
the  direction  of  the  rays  in  the  scene  beam,  produced  an 
effective  field  of  view  in  a  hologram  of  approximately  16 
degrees.   A  continuous  wave  helium-neon  laser  was  used  for 
alignment  of  the  O-switched  laser  and  the  optics.   A  Lauda 
constant  temperature  circulator  Model  N  controlled  by  an 
electronic  relay  type  R-10  was  used  in  conjunction  with  a 
Culligan  de-ionizer  to  maintain  the  laser  head  and  output 
etalon  at  a  constant  temperature  of  28.4  degrees  Centigrade. 
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In  Figure  h,    details  of  the  side  band  geometry  employed 
can  be  seen.   The  reference  beam  was  led  on  top  of  the  wind 
tunnel  by  means  of  the  two  mirrors  M  and  the  size  of  the 
beam  could  be  adjusted  by  means  of  lenses  U.   The  positions 
of  the  diffuser  plate  and  the  two  grids,  one  on  each  of  the 
wind  tunnel  side  walls,  are  also  shown  therein.   In  Figure 
5(a),  the  Korad  K-l  laser  and  the  smaller  alignment  laser 
are  shown  together  with  the  arrangement  of  the  optical 
equipment .   The  diffuser  plate  mounted  on  the  precision  X-Y 
slide  can  also  be  seen.   The  holographic  table  and  test  sec- 
tion was  completely  enclosed  in  a  wooden  box  to  enable  the 
experiments  to  be  conducted  in  the  daytime.   This  is  shown 
in  Figure  5(b).   Also  shown  therein  are  the  Korad  K-l  laser 
power  supply  and  shutter  control  assemblies,  the  Culligan 
de-ionizer  and  the  Lauda  constant  temperature  circulator. 

In  order  to  be  able  to  view  the  cone  at  angle  of  attack 
always  in  a  direction  perpendicular  to  its  axis  when  it  was 
rotated  in  the  wind  tunnel,  it  was  necessary  to  construct 
separate  grids  for  each  rotational  angle  at  which  the  holo- 
grams were  taken.   The  front  grid  used  for  the  82.5  degree 
rotational  position  of  the  cone  is  shown  in  position  in 
Figure  6.   The  grid  lines  can  be  seen  to  be  parallel  and 
perpendicular  to  the  cone  axis  for  this  particular  angular 
position.   The  orientation  of  the  grids  was  designed  to  be 
at  the  correct  angle  for  the  particular  rotational  position 
of  the  cone.   Since,  on  rotation,  the  axis  of  the  cone  was 
not  in  a  vertical  plane,  it  was  necessary  to  offset  the  two 
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grids  with  respect  to  each  other  in  such  a  manner  that  on 
lining  up  the  grids,  the  correct  viewing  position  perpendic- 
ular to  the  cone  axis  was  automatically  obtained.   The 
method  of  obtaining  the  orientations  and  the  relative  offset 
of  the  grids  is  shown  in  Figures  7(a)  and  7(b).   Since  the 
cone  spindle  was  bent  at  a  distance  of  0.7  inches  from  the 
base  of  the  cone,  the  various  rotational  positions  of  the 
cone  tip  are  indicated  in  the  end  view  in  Figure  7(a).   The 
corresponding  positions  of  the  cone  axis  in  the  side  and 
bottom  views  are  also  shown  therein.   The  orientation  of 
the  grid  lines  for  each  rotational  position  of  the  cone  as 
shown   in  Figure  7(b)  were  obtained  from  the  positions  of 
the  cone  axis  in  the  side  view  for  each  particular  case.   In 
order  to  be  normal  to  the  cone  axis,  lines  perpendicular  to 
the  orientation  of  the  cone  axis  in  the  bottom  view  were 
transcribed,  each  of  these  representing  a  line  of  sight. 
Since  the  distance  between  the  grids  was  7  7/8  inches,  the 
necessary  amount  by  which  the  two  grids  were  required  to 
be  offset  for  each  rotational  position  of  the  cone  tip  were 
obtained  as  shown  in  Figure  7(b). 

D.   THE  WIND  TUNNEL  MODEL 

The  cone  model  used  is  shown  in  Figure  8.   The  10  degree 
semi-apex  angle  cone  model  was  fabricated  out  of  stainless 
steel  and  had  pressure  ports  at  three  cross-sections,  the 
angular  position  of  the  ports  being  so  arranged  as  to  lie  on 
a  ray  from  the  apex  of  the  cone.   Since  a  fixed  mount  was 
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employed,  provision  was  made  to  attach  the  cone  to  spindles 
bent  to  various  angles  to  achieve  the  desired  angle  of  attack 
Stainless  steel  tubing  of  .036  inches  internal  diameter  was 
used  from  the  static  pressure  orifices  on  the  surface  of  the 
cone  up  to  the  base  of  the  model,  the  diameter  of  this  being 
adequate  to  provide  sufficient  time  response  of  the  Scani- 
valve.   Any  fifteen  out  of  the  thirty  ports  could  be  connec- 
ted at  a  time  to  the  Scanivalve  through  the  hollow  cone 
spindle  and  mount,  this  restriction  being  imposed  by  the 
internal  diameter  of  the  spindle. 

Rotation  of  the  model  about  the  horizontal  axis  of  the 
spindle  necessitated  drilling  4  holes  on  each  side  of  the 
wind  tunnel  side  walls  through  which  the  model  spindle  could 
be  unlocked  from  the  mount  and  rotated.   The  plexiglas  side 
walls  of  the  test  section  and  the  holes  therein  can  be  seen 
in  Figure  6.   The  collar  on  the  model  spindle  containing 
radial  holes  (also  visible  in  the  figure)  was  used  to 
rotate  the  cone  by  means  of  a  spike  arrangement  inserted 
through  one  of  the  holes  in  the  tunnel  side  walls.   The 
angle  of  rotation  was  established  by  lining  up  graduations 
on  the  spindle  sleeve  with  a  graduation  on  the  mount.   The 
model  and  spindle  used  for  the  axisymmetric  case  can  be  seen 
in  Figure  9.   The  pressure  ports  at  the  three  stations  and 
one  of  the  pressure  tubes  from  the  cone  passing  through  the 
hollow  spindle  are  visible  in  the  photograph.   After  the 
spindle  was  screwed  into  the  cone,  it  was  locked  in  position 
by  means  of  the  setscrew  shown. 
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III.   THEORETICAL  CONSIDERATIONS 

A.   EVALUATION  OF  THE  DENSITY  PROM  HOLOGRAPHIC  INTERFEROGRAMS 
1 .   The  Basic  Equation  of  Interferometry 

The  interferometry  principle  is  based  on  the  fact 
that  two  originally  coherent  rays  of  light  will  reinforce  or 
annul  each  other  at  a  screen  on  which  both  rays  are  pro- 
jected depending  on  their  relative  phase  there.   The  phase 
difference  is  directly  a  function  of  the  difference  in  the 
lengths  of  the  optical  paths  traversed  by  the  two  rays.   If 
the  optical  path  is  changed  by  an  amount  NX,  where  X  is  the 
wavelength  of  the  light  and  N  is  an  integer,  then  the  order 
of  interference  changes  by  an  amount  N,  that  is,  a  shift  of 
N  fringes  occurs.   The  change  in  optical  path  can  be  related 
to  the  density  through  the  index  of  refraction  and  the 

Gladstone-Dale  constant.   If  the  speed  of  light  in  a  medium 

c 

is  represented  by  —  where  c   is  the  speed  of  light  in 
^  J     n         o 

vacuum  and  n  is  the  index  of  refraction,  the  additional  tirre 

required  to  traverse  the  test  section  due  to  a  change  in  the 

index  of  refraction  from  n.,  to  n?  can  be  calculated  to  be 

equal  to  At  =  — (n?-  n, )  and  thus  the  change  in  optical  path 

o 
obtained  is  c  At,  that  is,  L(n2-n,).   The  fringe  shift,  g, 

which  is  equal  to  -r-  is  thus  given  by  L(rip-n,). 

For  a  given  substance  and  for  a  given  wavelength  of  light, 

the  index  of  refraction  is  a  function  of  density.   In  the 

case  of  gases,  since  the  speed  of  light  varies  only  slightly 
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from  that  in  vacuum,  it  can  be  represented  very  accurately 
by  the  first  two  terms  of  the  Taylor  series: 

n  =  /  +  p  L  (l) 

S 

where 

p   =  reference  density  at  0"C,  7^0  mm.  Hg . 

o 

3   =  dimensionless  constant  related  to  the  Gladstone-Dale 
constant  by: 

K-  £■ 

The  value  of  6  for  air  at  X    =  ?8°3  Angstroms  is  0. 0002^2, 
the  variation  with  wavelength  being  small.   The  fringe  shift 
relation  for  a  fixed  difference  in  index  of  refraction 
between  the  two  beams  becomes: 

*  '■ '  tC—f^O  (2> 

If  the  density  is  variable  in  the  test  section,  the  net 
change  in  optical  path  is  the  integrated  effect  along  the 
ray  which  gives  the  relationship: 


(3) 


where 


6J  -   *  £  (5) 
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ds  =  Incremental  distance  along"  a  ray 

z   =  a  plane  of  constant  z. 
c     ^ 

Equation  (3)  is  the  basic  equation  which  has  to  be 
inverted  to  obtain  the  unknown  density  from  the  known  fringe 
shift  values  obtained  from  an  Interferogram. 
2 .   The  Integral  Inversion 

The  integral  inversion  method  used  was  first  repor- 
ted by  C.  D.  Maldonado  et  al  in  1965  (1965;  Olsen,  1968)  and 
subsequently  used  by  R.  D.  Matulka  (Reference  3)  to  calculate 
the  density  variation  in  an  asymmetric  free  jet  from  the 
fringe  numbers  obtained  from  holographic  interferograms . 
The  procedure  Involves  the  representation  of  the  function 
f(x,y,z  )  by  a  complete  set  of  orthogonal  functions,  with 
the  unknown  coefficients  being  evaluated  by  use  of  the 
orthogonality  relationship  between  the  set  of  functions. 
The  set  of  functions  employed  have  a  region  of  orthogonality 
that  covers  the  entire  plane  and  have  the  property  that  they 
are  "invariant  in  form"  to  a  rotation  of  the  co-ordinate 
system.   The  co-ordinate  system  used  for  the  inversion  is 
shown  in  Figure  10  where  X  and  Y  represent  fixed  laboratory 
co-ordinates  and  X   and  Y   the  co-ordinates  in  which  the 
fringe  number  function  is  defined.   These  co-ordinates, 
therefore,  rotate  with  respect  to  the  fixed  co-ordinates  X 
and  Y  as  the  angle  of  view  through  the  test  section  is 
varied . 

In  operator  form,  equation  (3)  can  be  represented  as: 
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and  the  object  of  the  procedure  is  to  evaluate  f,  that  is, 
to  obtain: 

/<f*>*.  O   •  T"  y(5,%',  *c;  (7) 

where  T    represents  the  inverse  of  T. 

This  inversion  is  achieved  by  utilizing  a  pair  of 

orthogonal  polynomials  U~  ni  (ax, ay)  and  H  . _.  (ayT)  which  are 

m+2k  m+2k 

related  by  the  transform  relationship: 

T\U    (<**,  «v)  -c    C  a  _e /   .  /-/  ^yj  (8) 

where : 

rt,  ^K  C-  J  i-  «l         J  (  Q  ) 


and  H   are  the  Hermite  polynomials. 

The  function  f(x,y,z  )  is  first  expanded  in  a  series  of 

the  polynomials  U  . „.  as: 
^  J  m+2k 


(10) 


where  e   =  h   for  m  =  0  and  1  for  m  =  1,2,3, ,  **  is  an 

±m 
arbitrary  scale  factor,  and  C   ~,  are  the  unknown  expansion 

coefficients.   In  order  to  evaluate  these  coefficients,  an 

expression  for  the  function  g(5"?y^z  )  is  obtained  in  terms 

of  the  polynomials  by  transforming  f(x,y,z  )  in  accordance 

with  equation  (6).   Using  the  transform  relation  given  by 

equation  (8),  the  following  expression  is  obtained: 
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The  expansion  coefficients  can  then  be  evaluated  by 
applying  the  following  orthogonality  relationship  to 
Equation  ( 11 ) : 

/V""Vi",**$  /t.^  ^?»'J  «"*^/  - 


-/r 


that  is,  by  taking  the  scalar  product  of  equation  (11)  with 
.£     /7  ia     % )      and  using  equation  (12),  the  coefficients  of 
the  series  an  equation  (10)  are  obtained  as: 

Substitution  of  equation  (13)  into  equation  (10)  results 
■he  density  variation 


in  the  density  variation  being  given  by: 


.>,»*  **c  ^  (*»+**) i 


ReUJJ  1  ffa  I  O  *           " Zk  *      J      * r  ** 

+  m 

The  functions  IT  ?,  employed  are  complex  polynomials 

defined  as:                     .'/.,<-  ..  "7  '"Vx. 


U 

tr\   i-Zk. 


?d  as:  '/  r  7  '/ 


(15) 
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where  <J>  =  tan   J-r  and  L™  are  the  Associated  Laguerre 


T 

polynomials : 

6" 


(16) 


Substitution  of  equation  (15)  into  equation  (14)  yields: 

[b-£)  esc-*)  +  ^:^^^^^Y)(11) 

where : 


TT      oc 


B"    «>    ,  ff  m:*.**)  <•*(">*)%£>'>  4'<tf     (1 


8) 


Equations  (17),  (18),  and  (19)  are  the  basic  equations 
used  to  obtain  the  density  distribution  of  a  completely 
asymmetric  flow  field. 

3 .   The  Numerical  Procedure 

The  calculation  of  the  unknown  coefficients  P  , «,  (a) 

m+2k 

and  D   01  (a)  in  the  series  representation  of  f(x,y,z  )  as 
rrn~  c.  k.  c 

given  by  equation  (17)  is  not  analytically  possible  in  gen- 
eral since  the  function  g(yT,£,z  )  representing  the  variation 
of  fringe  number  is  an  experimentally  measured  quantity  and 
thus  not  explicitly  given  as  a  function  of  the  co-ordinates 
of  a  field  point.   If  the  analytical  procedure  is  to  be 
used  to  obtain  the  density  from  experimental  observations  of 
the  fringe  numbers,  therefore,  it  is  necessary  to  evaluate 
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the  double  integrals  in  equations  (18)  and  (19)  numerically. 
This  is  achieved  by  first  noting  that  for  every  g(y',£;,z  ) 
encountered  experimentally,  there  is  a  circular  domain  out- 
side which  the  fringe  distribution  is  zero  so  that  it  is 
possible  to  replace  the  limits  of  integration  of  +°°  and  -« 
in  equations  (18)  and  (19)  by  finite  values.   The  function 
g(y',£,z  )  is  then  approximated  by  constant  step  values  over 
elemental  areas  into  which  the  domain  is  partitioned.   This 
results  in  the  B  and  D  coefficients  being  represented  by  the 
double  series : 

S'V*)-/  7  ?($•  +  **;   *c-A*c)fecsO»S)<tsJ  H   C*x)  d%  (20) 

and 


I- 


T,+i  *■*' 


Using  the  derivative  formula  for  Hermite  polynomials, 
the  following  equations  result : 

{**•$.?  -^m#£.&h  -  'u.tr^J       (22) 

.[&*(»  f)  -  c^  (*,  f )  7/7/  («  *i  **)  -  H     C«  *;)  7    (23) 
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Since  it  is  not  possible  to  sum  over  an  infinite  number 
of  terms,  equation  (17)  is  replaced  by  the  sum  of  a  finite 
series : 


[Bml£) *>*<»>*)    +\&>  >-("<*)] eiJ""Y) 


(24) 


The  accuracy  of  the  representation  depends  on  the  values 
of  K3  M,  A£  and  Ax  chosen  and  convergence  of  the  series  is 
influenced  by  the  parameter  a. 

B.   THE  FLOW  AROUND  A  RIGHT  CIRCULAR  CONE 
1 .   Axisymmetric  Case 

The  problem  of  the  supersonic  flow  around  a  cone  has 
been  extensively  studied  by  several  investigators  ever  since 
it  was  first  opened  up  by  Busemann,  Bourquard  and  Taylor  and 
Maccoll.   The  basic  approach  to  the  problem  is  outlined  below 
but  a  detailed  description  may  be  found  in  References  12,  13, 
14  and  1?. 

Starting  with  the  equations  of  motion  of  inviscid 
flow  in  three  dimensions  expressed  in  spherical  co-ordinates: 


2 

.  o 


U  v}£      r     3£  *hi       r        ^-'    Jif    +  ,J_     ^J2     -t~     u-t--  used  C 


A^ 


(1) 
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and  the  equation  of  continuity 


L   (Kuf"v*j  +   i-  (f-*v  *u<-a)  +    4-  (f^^)   « 


J*. 


S  'J         7    ^°^^;  "-  °  (2) 


the  cylindrical  symmetry  of  the  cone  allows  one  to  eliminate 
the  velocity  component  w  and  all  terms  dependant  on  d>.      In 
addition,  application  of  the  condition  that  the  flow  be 
conical,  that  is,  no  variation  of  the  flow  parameters  with 
r,  reduces  equations  (1)  to: 

du.       _   ir    -    O 

gv  +  a,    +  .jL    <*£     t  o 

the  first  of  these  equations  asserting  that  the  flow  must  be 
irrotational .   The  equation  of  continuity  similarly  reduces 
to: 

Elimination  of  v  between  (3)  and  (4)  yields: 


Setting 


(6) 


Cf&  Uf        d£  ate 

the  following  ordinary  differential  equation  is  obtained: 

d\      +     u      ^        qY^-  *  tr-gcl  6.)  (?) 
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From  the  Bernoulli  equation: 


J       f 


JL  C 


(8) 


and  with  assumption  of  a  perfect  gas  satisfying  the  equation 
k 


r 


(9) 


the  velocity  of  sound  a  can  be  expressed  In  terms  of  the 
velocity  components  as: 

4  w- 


a 


do) 


The  density,  pressure  and  temperature  are  then  found  to 
vary  as  : 


/   v    P  /      *-  -  ^\   S-i 

(a)  L         ,    (c    -  u    -  v  ) 


r~, 


(11) 


(c).r    -f*/  •  c*- «'-*>) 


For  the  solution  of  the  governing  equation  (7)  to  be 
physically  meaningful,  it  is  necessary  to  satisfy  the  boun- 
dary condition  of  no  normal  velocity  at  the  cone  surface 
and  the  jump  conditions  at  the  shock  wave.   The  physical  con- 
ditions on  either  side  of  a  shock  wave  are  given  by  the  Ran- 
kine  Hugoniot  relations: 


2. 


2  * 


(&) 


a, 


m 


?  --" 
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J 


(12) 


where,  referring  to  Figure  11, 

a   =  Angle  between  the  normal  to  the  shock  wave  and  the 
streamlines  in  front  of  it. 

3   =  Angle  between  the  normal  to  the  shock  wave  and  the 
streamlines  immediately  behind  it. 

(J)   =  Angle  between  the  streamline  at  any  point  and  the 
radius  vector  from  the  vertex  of  the  cone. 

U   =  Freestream  velocity. 

a   =  Velocity  of  sound  in  the  freestream. 
x       / 

y    VFi 

—  =  p  /p 

z     ^w  ^o 

P-,  =  Pressure  of  the  undisturbed  gas  in  front  of  the 
shock  wave. 

p   =  Pressure  immediately  behind  the  shock  wave, 
w 

p   =  Stagnation  pressure  behind  the  shock  wave. 

The  boundary  condition  at  the  shock  wave  thus  reduces 
to  the  requirement  that  the  ratio  (x/z)  defined  by  the 
Rankine  Hugoniot  relations  be  equal  to  the  pressure  ratio 
(p/p  )  defined  by  equation  (11a).   This  together  with 
equations  (12)  defines  the  physical  characteristics  of  the 
shock  wave . 

Many  schemes  have  been  employed  in  order  to  numer- 
ically integrate  the  governing  differential  equation  and 
to  satisfy  the  boundary  conditions  to  yield  the  complete 
solution  of  the  flow  field.   Results  of  some  of  these  are 
tabulated  in  References  6  and  14.   For  the  purposes  of  this 
report,  the  more  recent  results  in  Reference  6  have  been 
used . 
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2 .   Flow  Around  Cones  with  Small  Yaw 

The  governing  differential  equation  for  the  flow 
around  a  yawed  cone  is  developed  on  similar  lines  as  that 
for  axisymmetric  flow  with  the  difference  that  in  this  case, 
the  velocity  components,  pressure  and  density  are  expressed 
as  a  series  representation  about  the  axial  flow  values.   A 
detailed  description  of  the  procedure  may  be  found  in 
References  13,  15  and  17.   An  examination  of  the  results 
reveals  that  although  the  flow  is  no  longer  irrotational , 
the  flow  field  can  still  be  regarded  as  conical.   Further- 
more, the  shock  wave  in  this  case  is  still  a  cone  of  the 
same  semi-apex  angle  as  in  the  no-yaw  case  with  the  dif- 
ference that  the  shock  wave  is  itself  yawed  at  an  angle 
which  is  proportional  to  the  angle  of  yaw  of  the  cone  with 
respect  to  the  freestream. 

In  Reference  6,  the  method  employed  to  obtain  the  flow 
properties  is  basically  one  of  Iteration.   An  estimate  is 
made  of  the  shape  of  the  shock  wave  and  this  shape  is 
improved  in  such  a  way  that  the  normal  velocity  at  the  sur- 
face of  the  cone  is  made  close  to  zero.   The  results  are 
presented  in  the  co-ordinate  system  shown  in  Figure  12. 
Values  of  u,  v,  w,  — ,  and  p  at  various  azimuthal  angles  cj> 
are  given  for  various  non-dimensional  distances  between  the 
body  and  the  chock  wave.   Such  results  are  tabulated  for 
cones  of  different  half-cone  angles,  Mach  number  and  angle 
of  attack  at  which  the  flow  can  no  longer  be  considered 
conical.   For  the  case  of  a  10  degree  half-angle  cone, 
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this  limiting  angle  of  attack  is  10  degrees  at  a  Mach 
number  of  3.0  and  11  degrees  at  a  Mach  number  of  2.0. 


26 


IV.   DESCRIPTION  OF  THE  EXPERIMENT 

A.   NUMERICAL  SIMULATION  OF  THE  EXPERIMENT 

The  suitability  of  the  computer  program  HOLOFER  for 
handling  cases  with  discontinuities  in  slope  had  been  inves- 
tigated by  means  of  various  test  function  inputs  into  the 
program  In  Reference  3.   The  accuracy  of  the  inversion  In 
these  cases  was  reported  to  be  within  3.8  percent.   In  order 
to  verify  the  applicability  of  the  program  to  the  present 
investigation,  known  values  of  the  density  along  one  line 
of  sight,  obtained  from  Reference  6  at  a  particular  cross- 
section  in  the  flow,  were  fed  in  as  input  and  Mode  1 
operation  of  the  computer  program  utilized.   In  this  mode, 
the  program  calculates  the  fringe  number  array  that  would 
result  from  the  specified  density  distribution  and  then 
uses  this  data  as  the  input  for  inversion  to  obtain  the 
density  distribution.   One  of  the  primary  objectives  of  this 
investigation  was  to  determine  the  effect  of  opaque  objects 
introduced  into  the  flow  field  and  the  presence  of  the  shock 
wave  on  the  resulting  inversion.   It  was  found  that  the 
effect  of  the  presence  of  the  cone  could  be  reduced  by  the 
introduction  of  a  fictitious  density  distribution  in  the 
area  occupied  by  it  of  a  constant  value  equal  to  the  density 
at  the  cone  surface. 
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B.   EXPERIMENTAL  PROCEDURE 

1 .   Laboratory  Techniques 

It  has  been  shown  in  Reference  10  that  the  trans- 
mitted intensity  from  a  diffuse  glass  falls  off  below  30 
percent  of  the  incident  intensity  beyond  viewing  angles  of 
±8  degrees.   Thus  this  represents  a  limit  on  the  field  of 
view  obtainable  from  a  single  hologram.   For  the  flow  around 
a  right  circular  cone  at  Incidence,  there  is  planar  symmetry 
and  hence  a  90  degree  field  of  view  is  required  for  complete 
coverage  of  the  total  flow  field.   This  coverage  can  be 
obtained  by  the  simultaneous  production  of  a  number  of 
holograms  in  various  azimuthal  positions  as  employed  in 
Reference  3  or  by  a  rotation  of  the  viewing  angle  between 
each  production  of  a  single  hologram.   Taking  advantage  of 
the  steadiness  of  the  flow  fields  investigated,  holograms 
covering  various  sectors  of  the  flow  field  were  obtained 
by  rotating  the  cone  to  six  positions  in  a  90  degree  sector, 
one  hologram  being  obtained  in  each  case. 

Complete  angular  coverage  of  the  flow  field  was 
obtained  by  a  rotation  of  the  cone  and  spindle  arrangement 
in  the  mount  by  an  angle  of  15  degrees  between  each  hologram 
Since  it  was  desired  to  obtain  the  flow  field  at  the  same 
cross-section  for  the  various  angles  of  rotation,  it  was 
necessary  to  obtain  the  holographic  interf erograms  in  such 
a  way  that  the  line  of  sight  was  always  perpendicular  to 
the  axis  of  the  cone.   For  this  purpose,  a  pair  of  grids 
were  used,  one  on  each  wall  of  the  test  section,  in  which 
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the  lines  were  designed  to  be  parallel  and  perpendicular 
to  the  cone  axis  for  a  particular  rotational  position  of 
the  cone.   The  two  grids  were  offset  laterally  with  respect 
to  each  other  in  order  to  account  for  the  rotation  of  the 
cone  axis  away  from  the  central  plane  when  the  cone  was 
rotated.   On  production  of  the  holograms,  the  variable 
viewing  property  of  light  field  holograms  was  utilized  to 
line  up  the  two  grids  thereby  ensuring  that  the  line  of 
sight  was  normal  to  the  axis  of  the  cone.   This  would  have 
involved  a  considerable  amount  of  complexity  and  effort  if 
a  conventional  Mach-Zehnder  interferometer  had  been  employed 
for  the  same  purpose  . 

Since  the  coherence  length  of  the  pulsed  laser  was 
only  about  2  cms.  at  the  output  used,  it  was  necessary  to 
maintain  the  optical  path  lengths  in  the  reference  and  scene 
beams  nearly  equal.   Since  the  scene  beam  path  traversed  4 
inches  of  plexiglas  while  the  reference  beam  passed  through 
two  thick  plano-convex  lenses,  it  was  necessary  to  compen- 
sate for  this  by  making  the  reference  beam  path  3-5  inches 
longer  than  that  of  the  scene  beam.   A  reference  to  scene 
beam  ratio  of  4:1  was  found  to  yield  good  holograms. 

The  holograms  were  obtained  on  Agfa-Gevaert  8E-75 
holographic  plates,  V  x  V  in  size.   Development  was  for  5 
minutes  in  Kodak  D-19  developer,  thirty  seconds  in  an  acetic 
acid  stop  bath  of  standard  dilution,  five  minutes  in  standard 
fixer,  one  minute  water  wash  followed  by  immersion  in  a 
Kodak  PhotoPlo   wetting  agent  prior  to  drying.   Reconstruction 
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of  the  scene  was  achieved  by  illuminating  the  hologram  with 
the  light  from  a  helium-neon  laser  operating  at  a  wavelength 
of  6328  Angstroms. 

2 .   Photographic  Techniques 

To  obtain  reconstruction  of  the  image  from  the  holo- 
gram, various  methods  can  be  employed.   In  the  usual  method 
where  the  hologram  is  re-illuminated  by  a  beam  similar  to 
the  original  reference  beam  as  shown  in  Figure  13,  each 
point  on  the  photograph  is  produced  by  a  series  of  non- 
parallel  rays  that  appear  to  originate  from  various  sources 
on  the  diffuse  glass  that  was  used  to  construct  the  hologram, 
An  almost  parallel  set  of  rays  can  be  selected  by  the 
positioning  of  a  small  aperture  at  the  focal  plane  of  the 
imaging  lens  as  shown  in  Figure  14.   A  similar  effect  is 
achieved  by  illuminating  the  hologram  by  a  conjugate  beam 
of  small  diameter  as  depicted  in  Figure  15.   Because  of  the 
small  diameter  of  the  reconstruction  beam,  the  illuminated 
portion  of  the  hologram  represents  a  small  aperture  and  thus 
the  image  produced  has  a  large  depth  of  field.   This  was  of 
considerable  advantage  in  the  production  of  the  interfero- 
grams  from  the  holograms  since  it  enabled  the  front  and  rear 
grids,  the  cone  and  the  fringes  to  be  simultaneously  pro- 
jected on  a  screen  and  thus  made  the  task  of  obtaining  the 
correct  viewing  angle  and  the  interferogram  quite  easy  to 
achieve.   A  photograph  of  the  scene  was  obtained  by  focusing 
the  camera  at  the  desired  plane  in  the  test  section,  usually 
in  such  a  manner  so  as  to  be  in  the  plane  of  the  fringes. 
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The  lines  of  sight  recorded  on  the  photograph  in  this  case 
represent  the  almost  parallel  pencil  of  rays  from  the 
diffuser  that  pass  through  the  effective  aperture  in  the 
hologram. 

3 .   Data  Reduction 

To  obtain  photographic  interferograms ,  the  hologram 
was  illuminated  by  the  conjugate  reference  beam  using  a 
continuous  wave  laser  operating  at  a  wavelength  of  6328 
Angstroms.   A  camera  back  with  viewing  screen  was  placed  in 
the  position  of  the  real  image  as  shown  in  Figure  13.   By 
lining  up  the  appropriate  horizontal  and  vertical  lines  of 
the  front  and  rear  grids,  the  viewing  angle  corresponding 
to  the  particular  line  of  sight  desired  was  achieved.   As 
the  spacing  of  the  lines  on  each  grid  was  \   inch  and  the 
two  grids  were  physically  7  7/8  inches  apart,  a  viewing 
angle  of  3° 38'  was  obtained  by  lining  up  a  line  on  one  of 
the  grids  with  the  adjacent  one  on  the  other.   A  f 7 . 7  lens 
with  full  aperture  and  an  exposure  time  of  4  seconds  was 
used  to  obtain  workable  holograms  on  Poloroid  Type  55  P/N 
film.   The  reduction  of  fringe  shift  was  obtained  by  project- 
ing the  negative  in  a  photoenlarger  and  by  tracing  the 
fringes,  the  cone  surface,  the  shock  wave  and  the  grid  line 
at  the  desired  station  at  which  the  reduction  was  to  be 
performed,  on  a  sheet  of  paper.   Some  judgement  was  usually 
necessary  in  locating  the  centre  of  the  fringes  but  in 
general  the  light  fringes  were  found  to  be  narrower  than  the 
dark  fringes  and  a  location  of  the  fringes  to  within  ±0.1  of 
the  fringe  spacing  in  the  freestream  could  be  obtained.   The 
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fringes  in  the  region  outside  the  shock  were  extended  towards 
the  cone  surface  and  the  fringe  shift  read  at  the  points  of 
intersection  of  the  displaced  fringes  with  a  number  of  rays 
from  the  cone  vertex.   The  assumption  of  locally  conical 
flow  in  the  region  within  two  to  three  fringes  on  either  side 
of  the  section  was  made  in  order  to  obtain  the  fringe  dis- 
placements at  the  section  from  those  at  the  points  of 
intersection  of  the  rays  from  the  cone  vertex  with  the 
fringes  on  either  side  of  the  section.   The  distance  between 
the  cone  tip  and  the  section  at  which  reduction  was  per- 
formed was  measured  from  the  projected  image  and,  since  the 
actual  distance  between  these  points  was  known,  it  was 
possible  to  calculate  the  factor  by  which  the  projected 
image  was  magnified  as  compared  with  its  true  size.   This 
magnification  factor  was  then  used  to  determine  the  actual 
radial  distances  from  the  cone  axis  at  which  the  fringe 
shifts  were  obtained.   The  radius  of  the  inversion  circle 
was  selected  such  that  the  maximum  distance  of  the  shock 
from  the  cone  axis  was  at  95  percent  of  the  radius  of  the 
inversion  circle.   Prom  the  data  so  obtained,  the  radial 
variation  of  fringe  number  was  plotted  and  a  smooth  curve 
drawn  through  the  points.   The  fringe  number  at  201  equidis- 
tant points  was  then  read  off  from  the  curve  and  utilized 
as  the  input  into  the  computer  program  HOLOPER  in  Mode  3- 
This  mode  utilizes  raw  data  taken  at  the  proper  interval 
and  directly  read  into  the  G  array  in  the  program  by  Sub- 
routine READ.   Further  details  about  the  use  of  this  computer 
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program  are  outlined  in  Appendix  C.   For  the  axi symmetric 
case,  it  was  necessary  to  feed  in  the  fringe  data  along  one 
line  of  sight  only  and  to  obtain  the  density  variation  along 
one  radial  line.   For  the  asymmetric  case,  fringe  data  was 
fed  in  along  six  lines  of  sight  in  a  90  degree  field  of  view 
and  the  inverted  density  field  obtained  along  9  radial  lines 
spanning  a  180  degree  field  of  view  on  one  side  of  the  axis 
of  planar  symmetry  of  the  flow.   A  typical  reduction  of  an 
interf erogram  to  obtain  the  radial  variation  of  fringe 
number  at  a  section  is  shown  in  Appendix  A. 
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V.   EXPERIMENTAL  RESULTS  AND  DISCUSSION 

A.   DENSITY  FIELD  OBTAINED  EXPERIMENTALLY 
1 .   Computer  Simulation  Results 

The  density  distribution  obtained  from  the  AGARD  167 
tables  for  the  axisymmetric  case  was  used  as  a  test  function 
in  Computer  program  HOLOPER .   The  fringe  distribution  cal- 
culated from  this  density  distribution  was  inverted  to  obtain 
the  original  density  distribution.   A  plot  of  the  original 
density  distribution  and  that  inverted  by  the  program  is 
shown  in  Figure  16.   Since  the  Mach  number  for  this  inversion 
was  2.84,  It  was  necessary  to  first  interpolate  the  13  values 
obtained  from  Reference  6  for  a  Mach  number  of  2.84,  and 
then  to  obtain  from  these,  the  values  at  201  equidistant 
points  along  a  line  of  sight  within  the  region  of  the 
inversion  circle.   This  was  done  using  Computer  Program  1. 
As  a  side  benefit  obtained  from  a  use  of  this  program,  the 
discontinuities  in  slope  in  the  density  distribution  at  the 
cone  surface  that  would  otherwise  have  resulted  by  assuming 
a  constant  value  In  the  region  occupied  by  the  cone,  were 
eliminated.   Various  values  of  K,  M  and  interval  size  were 
experimented  with  to  obtain  the  minimum  deviation  of  the 
inverted  density  distribution  from  the  input  value  in  the 
region  close  to  the  shock  wave.   Increasing  the  value  of  K 
above  500  was  found  to  have  little  effect  on  the  accuracy, 
the  maximum  variation  of  the  inverted  density  being  1.5 
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percent  of  the  input  density  value  at  the  point  of  inversion. 
The  value  of  a  was  found  to  affect  the  total  number  of  terms 
in  the  series  representation  of  the  B  and  D  coefficients 
required  for  convergence  within  the  value  of  e  specified.   A 
value  of  K  =  2000,  M  =  5,  a  =  2.5  and  an  interval  size  of 
0.01  was  found  to  yield  the  minimum  effect  of  the  shock 
discontinuity . 

2 .   Experimental  Results 
a.   Axisymmetric  Case 

A  photograph  of  the  holographic  interferogram 
obtained  is  shown  in  Figure  17(a).   The  fringe  data  obtained 
from  a  reduction  of  the  two  sides  of  the  cone  is  shown  in 
Figure  17(b).   Since  the  flow  was  axisymmetric,  an  average 
of  the  two  readings  was  taken  and  a  curve  plotted.   The 
values  at  201  equidistant  points  was  read  off,  this  being 
input  into  Subroutine  READ  of  computer  program  HOLOFER  as 
the  G  array.   The  density  distribution  obtained  from  the 
output  of  the  inversion  is  shown  in  Figure  18.   Also  shown 
plotted  is  the  density  field  obtained  from  the  AGARD  167 
tables  interpolated  for  a  Mach  number  of  2.84  by  computer 
program  2,  this  being  the  Mach  number  of  the  tunnel  obtained 
from  pressure  measurements.   The  pressure  data  obtained 
from  the  viscicorder  is  shown  in  Figure  19.   The  value  of 
Mach  number  calculated  from  these  pressure  measurements  com- 
pares favorably  with  the  value  obtained  in  earlier  experi- 
ments from  shadowgraphs  of  the  flow  in  the  tunnel.   A 
comparison  of  the  density  distributions  from  the  AGARD  tables 
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and  those  obtained  from  the  experiment  show  good  agreement , 
the  maximum  variation  in  overall  density  being  8  percent  of 
the  theoretical  density  value, 
b.   Asymmetric  Case 

The  views  along  which  holograms  were  taken  for 
the  cone  at  10  degrees  angle  of  attack  are  shown  in  Figure 
20  and  photographs  of  the  interferograms  obtained  are  shown 
in  Figures  21  to  26.   Except  for  the  67.5  degree  and  82.5 
degree  angles  of  view,  the  perturbation  to  the  flow  on  the 
leeward  side  of  the  cone  was  insufficient  to  enable  either 
the  shock  position  or  the  fringe  displacements  on  this 
portion  of  the  flow  field  to  be  obtained.   The  fringe  shifts 
in  this  region  were,  therefore,  taken  as  being  zero.   A 
comparison  of  the  shock  wave  position  obtained  experimental- 
ly from  the  interferograms  with  that  from  the  AGARD  167 
tables  interpolated  for  a  Mach  number  of  2.84  is  shown  in 
Figure  27.   The  fringe  data  were  input  into  computer  program 
HOLOFER  starting  with  the  82.5  degree  fringe  distribution 
in  order  to  be  consistent  with  the  plane  of  symmetry  for 
which  that  program  was  written.   The  resulting  density 
distribution  output  from  the  program  was  along  the  lines  of 
sight  indicated  in  Figure  28.   A  comparison  of  the  experi- 
mentally obtained  radial  distributions  with  those  from 
the  AGARD  167  tables  interpolated  for  a  Mach  number  of  2.84 
are  shown  in  Figures  29  to  35. 
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B.   ERROR  ANALYSIS 

The  errors  in  the  final  solution  are  mainly  due  to  errors 
in  the  fringe  data  input  into  the  computer  program.   One  of 
the  major  sources  of  error  was  due  to  the  difficulty  in 
obtaining  the  slope  of  the  fringe  lines  in  the  freestream  to 
a  very  high  degree  of  accuracy.   Prom  the  results  of  the 
experiments  conducted  on  a  10  degree  half  angle  cone  In 
Reference  18,  it  was  apparent  that  the  flow  in  regions  close 
to  the  tip  of  the  cone  would  be  affected  by  the  boundary 
layer  on  the  cone  whereas  sections  near  the  base  of  the  cone 
would  be  affected  by  the  wake  and  expansion  of  the  flow.   In 
order  to  obtain  a  flow  which  was  locally  nearly  conical,  it 
was  thus  necessary  to  reduce  the  interf erograms  at  a  section 
approximately  60  to  70  percent  of  the  cone  length  downstream 
of  the  tip.   At  this  section,  however,  the  shock  wave  was 
about  1  inch  from  the  axis  of  the  cone  and  considering  that 
the  fringes  in  a  region  about  \    inch  adjacent  to  the  tunnel 
walls  were  affected  by  the  boundary  layer,  this  left  a  region 
of  about  3/^  inch  in  which  the  fringes  of  the  flow  in  the 
freestream  were  present.   It  was  thus  difficult  to  obtain 
the  slope  of  the  fringes  very  accurately  and  it  is  estimated 
that  errors  of  approximately  ±1  mm.  in  the  fringe  data  along 
any  fringe  could  have  resulted.   However,  since  the  fringe 
data  along  the  section  was  obtained  by  interpolation  between 
different  fringes  at  varying  radial  distances  from  the  cone 
axis,  this  would  have  resulted  in  a  smoothing  of  the  errors 
to  about  ±0.5  mm.   Since  the  maximum  fringe  shifts  obtained 
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were  in  the  region  of  8  mm.,  an  error  of  ±0.5  mm.  corresponds 
to  an  error  of  approximately  ±8  percent  in  the  magnitudes  of 
the  fringe  shifts  obtained. 

The  errors  between  the  experimental  and  theoretical 
density  curves  in  Figure  18  reflect  the  above  estimates  of 
error  quite  well.   The  experimentally  obtained  density 
distribution  is  generally  about  5%   higher  than  the  theoreti- 
cal curve  and  this  rises  to  a  maximum  of  8%    close  to  the 
cone  surface.   Similar  errors  are  apparent  in  general  in 
Figures  29  to  35.   However,  in  these  curves  the  maximum  error 
is  somewhat  higher  probably  due  to  the  fact  that  reduction 
was  performed  at  a  section  closer  to  the  base  of  the  cone 
than  for  the  axisymmetric  case.   In  so  far  as  the  general 
magnitude  of  the  density  is  concerned,  the  experimentally 
obtained  distribution  follows  the  theoretical  one  fairly 
well.   Except  for  the  density  distributions  for  cj>  =  0°  and 
cf)  =  22.5°,  the  density  starts  to  rise  above  the  freestream 
value  at  a  normalized  radial  distance  from  the  cone  axis 
of  about  0.7,  then  assumes  a  nearly  constant  value  before 
rising  steeply  again  at  a  radial  distance  closer  to  the 
cone  axis.   This  appears  to  be  inconsistent  with  the  data 
Input  into  the  program  in  which  the  radial  distance  of  the 
shock  from  the  cone  axis  is  seen  to  decrease  with  increasing 
values  of  (J). 
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V.       SUMMARY 

A.  CONCLUSIONS 

The  finite  fringe  method  for  the  production  of  holo- 
graphic interferograms  has  been  applied  successfully  to 
the  determination  of  the  three-dimensional  density  distri- 
bution of  the  flow  around  a  10  degree  half  angle  cone.   The 
method  has  been  demonstrated  to  yield  results  within  an 
accuracy  of  8  percent  for  the  axisymmetric  case  and  within 
13  percent  for  the  asymmetric  case.   The  magnitude  of  the 
errors  are  mainly  a  result  of  the  constraints  imposed  by 
the  existing  wind  tunnel  size  and  the  fact  that  the  cone 
model  used  produced  perturbations  in  the  flow  field  of  a 
fairly  small  magnitude  so  that  small  errors  in  measurement 
were  reflected  as  relatively  large  percentage  errors.   The 
computer  program  HOLOFEP  has  been  evaluated  for  the  effect 
of  the  shock  and  the  cone  and  found  to  yield  good  results. 
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FIGURE  2.    WIND  TUNNEL    TEST    SECTION  AND    PRESSURE 
MEASURING    EQUIPMENT 
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FIGURE  3.     PRESSURE    RECORDING   EQUIPMENT,  CONTROL 
BOX,  AND  POWER  SUPPLIES 
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FIGURE  5a.    HOLOGRAPHIC    TABLE    LASER    CONTROL    BOX 
AND    TEMPERATURE   CONTROL   EQUIPMENT 
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FIGURE  5b.     DETAILS  OF  THE    HOLOGRAPHIC  ARRANGEMENT 
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FIGURE  6.    VIEW  OF  TEST    SECTION     SHOWING    DETAILS   OK 
MODEL,  MOUNT,  AND   GRID 
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Orientation  of  Grid  Lines  for  Various  Positions 
of  the  Cone  Top 
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FIGURE  7b.  ORIENTATION  OF  GRID  LINES  AND  OFFSET  OF 
THE  GRIDS 
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FIGURE  9.    THE  CONE  AND     SPINDLE    ARRANGEMENT 
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FIGURE   10.    CO-ORDINATE      SYSTEM    USED    FOR   THE     INVERSION 
OF  FRINGE    NUMBER    TO    DENSITY 
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Streamlines 


FIGURE  II.    CO-ORDINATE    SYSTEM    USED  FOR  THE    GAS 
DYNAMIC  EQUATION 
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FIGURE  13.    EFFECT  OF  APERTURE    SIZE    FOCUS  PLANE 

POSITION   ON  THE  PENCIL   SIZE  OF  RAYS  ABOUT 
A   LINE  OF  SIGHT   RECORDED    BY  CAMERA 
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FIGURE   14.     SPATIAL     FILTERING    TECHNIQUE    FOR    SELECTING 
PHOTOGRAPH    OF  CONSTANT  ANGLE  LINES  OF 
LIGHT 
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FIGURE  15.    LENSLESS    PHOTOGRAPHIC  TECHNIQUE    USING  A 

CONJUGATE   REFERENCE  BEAM   OF  SMALL  DIAMETER 
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FIGURE  17a.  HOLOGRAPHIC    INTERFEROGRAM    OBTAINED    FOR 
THE  AXI- SYMMETRIC   CASE 
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FIGURE  17b.    RADIAL    FRINGE    DISTRIBUTION   FROM  A  REDUCTION   OF 
THE    HOLOGRAPHIC    INTERFEROGRAM    FOR    THE 
AX  I-  SYMMETRIC    CASE 

60 


3 

cvi 
ii 


o 


xwwwwv 


\\\\\\\\\\\VA 


O) 


CO 


x 
< 

c 
o 
o 

E 
o 


0> 

o 

CO 


o 


o 

c 
o 

«:  I 


c 
o 


ro 


LU 
CO 

5  < 

O  <-> 

cr 
Q  F 

UJ  UJ 

?  £ 

O  J. 


I- 

3 
CD 


CO 
Q 

CO 


O 

I- 

m 

tr 

h- 
co 

Q 

CO 
UJ 
Q  Q 

<  CD 

5  < 

o  * 


2?  |  " 

K!  <  £ 

^  q:  lu 

o  2  cl 

iT  °  X 

U-  O  UJ 


ro 


CM 


<N 


00 


(     d/d)    OIJDy    A|jSU9Q 


61 


Pressure  Ports 
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FIGURE  19.     PRESSURE     TRACE     OBTAINED    FROM     THE     VISCICORDER 
FOR    THE     AXI -SYMMETRIC     CASE. 
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Viewing  Positions 
82.5° 


22.5 


7.5' 


Views  Obtained 
in  Holograms 


22.5 


7.5* 


Cone  at  10°  Angle    of 
Attack  (Tip  Below  the 
Axis) 


Axis  of  Planar   Symmetry 


FIGURE  20.   VIEWS    ALONG    WHICH    HOLOGRAMS    WERE    OBTAINED 
FOR    THE    ASYMMETRIC    CASE 


63 


FIGURE   21.     HOLOGRAPHIC    INTERFEROGRAM    FOR    7.5° 
ANGLE  OF  VIEW 


6h 


FIGURE   22.  HOLOGRAPHIC     INTERFEROGRAM    FOR    22.5 
ANGLE  OF  VIEW 
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FIGURE    23.   HOLOGRAPHIC     INTERFEROGRAM    FOR      37.5° 
ANGLE    OF  VIEW 
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FIGURE   24.    HOLOGRAPHIC      INTERFEROGRAM      FOR    52.5° 
ANGLE  OF  VIEW 
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FIGURE   25.     HOLOGRAPHIC     INTERFEROGRAM    FOR    67.5° 
ANGLE  OF   VIEW 
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FIGURE  26.    HOLOGRAPHIC    INTERFEROGRAM    FOR    82.5° 
ANGLE   OF   VIEW 
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Cone 


tf  Experiment 
o  From  AGARD 


7.5°-  15°    Experimental 
0°-22.5°    AGARD 


FIGURE   27.    COMPARISON    OF   THE   POSITION  OF  THE    SHOCK    WAVE 
OBTAINED    EXPERIMENTALLY    WITH   THAT   FROM    AGARD 
167    FOR   M  =  2.84  ,  ASYMMETRIC   CASE 
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Axis  of  Planar   Symmetry 


180°     ,5Z5° 


FIGURE  28. 


LINES  OF  SIGHT  ALONG  WHICH  THE  DENSITY 
FIELD  WAS  OBTAINED  AFTER  INVERSION,  ASYM- 
METRIC   CASE 
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TABLE  I 

VARIATION  OP  FRINGE  SHIFT  WImH  DISTANCE 
ALONG  RAYS  FROM  THE  CONE  VERTEX 

Fringe  shift  (mm.) 

Distance  from  cone  vertex  along  a  ray  (ins.) 


Ray 

Fringe 
Ta' 

Fringe 

Fringe 
'c' 

Fringe 

Fringe 
'eT 

Section 

X 

y 

X 

y 

X 

y 

X 

y 

X 

y 

X 

y 

1 

0.0 

4.50 

0.0 

4.63 

0.0 

4.75 

0  .0 

4.89 

0.0 

5.0 

4.48 

2 

2.0 

4.37 

2.1 

4.40 

1.9 

4.63 

1.7 

4.77 

1.0 

4.92 

4.44 

3 

3.3 

4.26 

5.5 

4.40 

3.6 

4.50 

3.2 

4.64 

2.1 

4.78 

4.4n 

4 

4.1 

4.17 

4.4 

4.29 

4.4 

4.41 

4.1 

4.53 

3.5 

4.68 

4.36 

5 

4.9 

4.08 

4.9 

4.21 

5.4 

4.32 

5.0 

4.43 

4.7 

4.58 

4.32 

6 

5.5 

4.02 

5.6 

4.13 

6.1 

4.24 

5.5 

4.35 

5.6 

4.48 

4.28 

7 

6.0 

3-95 

6.1 

4.08 

6.6 

4.18 

6.0 

4.29 

6.4 

414 

4.26 

8 

6.4 

3.91 

6.8 

4.03 

7.2 

4.13 

6.7 

4.22 

6.9 

4.34 

4.23 

9 

6.6 

3.87 

7.2 

3.97 

7.7 

4.08 

6.6 

4.18 

7.4 

4.28 

4.21 

10 

6.6 

3.84 

7.3 

3.92 

7.8 

4.03 

6.6 

4.14 

7.7 

4.23 

4.18 

11 

6.5 

3.80 

7.1 

3-91 

7.2 

4.0 

6.5 

4.10 

7-7 

4.19 

4.15 

12 

6.3 

3.79 

6.4 

3-90 

7.1 

3.99 

6.3 

4.08 

7.3 

4.17 

4.i4 
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Ray 


TABLE  II 

RADIAL  DISTANCE  FROM  THE  CONE  AXIS  OP  THE 
INTERSECTION  OF  RAYS  VI^H  THE  SECTION 


Distance  from        Normalized  distance 
cone  axis  from  cone  axis 


1  1.83  0.95 

2  1.72  O.89 

3  1.63  0.846 

4  1.53  0.794 

5  1.43  0.742 

6  1.33  0.69 

7  1.23  O.639 

8  1.12  0.581 

9  1.02  0.53 

10  0.92  0.478 

11  0.82  0.426 

12  0.72  0.374 
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TABLE  III 

CALCULATION  OP  THE  RADIAL  VARIATION 
OF  FRINGE  NUMBER  AT  THE  SECTION 

Average  fringe  spacing  in  the  freestream  =  2.775  nrn 


Ray 


Normalized  distance 
from  cone  axis 


Fringe  shift 
( mm .  ) 

0. 

n 

2. 

n8 

3. 

5 

4. 

,42 

5- 

,4 

5. 

•  92 

6, 

,16 

6, 

.69 

6 

.73 

6. 

.87 

7 

.22 

6 

.80 

Fringe  number 


1 

0.95 

2 

0.89 

3 

0.846 

4 

0.794 

5 

0.742 

6 

0.69 

7 

0.639 

8 

0.581 

9 

0.53 

10 

0.478 

11 

0.426 

12 

0.374 

0. 

0 

0. 

75 

1. 

,26 

1, 

■  59 

1, 

.Q46 

2 

.13 

2 

.22 

2 

.41 

2 

.425 

2 

.476 

2 

.60 

2 

.45 
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APPENDIX  A 
REDUCTION  OF  AN  INTFRFEFOGFA^  TO  OPTAIN  FRINGE  SHIFT  DATA 

The  reduction  of  only  one  side  of  the  interferogram  ob- 
tained for  the  axi symmetric  case  is  indicated  as  an  illus- 
tration of  the  procedure  employed.   .After  projection  of  the 
negative,  the  cone  surface,  the  shock  waves,  the  prid  line 
at  the  section  concerned  and  the  fringes  from  the  projected 
image  were  traced  on  a  sheet  of  paper  as  shown  in  Figure  36. 
The  distance  between  the  cone  surface  and  the  shock  at  the 
section  was  divided  into  a  convenient  number  of  parts  and 
rays  drawn  from  the  tip  of  the  cone  through  these  points. 
The  fringe  shift  distances  at  the  point  of  intersection  of 
each  displaced  fringe  and  the  rays  were  measured  by  means  of 
a  7X  PEAK  scale  magnifier  as  well  as  the  average  fringe 
spacing  in  the  freestream.   ^he  radial  distance  from  the 
cone  axis  to  the  intersection  of  each  of  the  rays  with  the 
section  was  also  measured  as  also  the  distance  along  each 
ray  to  the  intersection  with  each  fringe  and  with  the 
section  concerned.   A  tabulation  of  these  results  is  shown 
in  Tables  1,  2,  and  3.   For  each  ray,  the  fringe  shift  dis- 
tance was  then  plotted  against  the  distance  along  the  ray 
as  in  Figure  37,  and  from  these  curves  the  fringe  shift 
distance  at  the  section  was  obtained.   These  fringe  shift 
distances  were  divided  by  the  average  fringe  spacing  in  the 
freestream  to  obtain  the  fringe  numbers  at  the  various 


points  concerned.   Prom  a  measurement  of  the  distance  from 
the  tip  of  the  cone  to  the  section  on  the  projected  image 
and  comparison  with  the  known  distance,  the  magnification 
of  the  image  with  respect  to  the  actual  conditions  existing 
in  the  wind  tunnel  was  determined.   The  radius  of  the  inver- 
sion circle  was  obtained  by  assuming  the  distance  of  the 
shock  from  the  cone  axis  to  be  at  95%    of  the  radius  of  the 
inversion  circle.   The  distance  from  the  cone  axis  of  the 
various  intersections  of  the  rays  from  the  cone  vertex  and 
the  section  were  obtained  as  a  fraction  of  the  inversion 
circle  radius.   The  fringe  number  at  these  points  was  then 
plotted  against  the  non-dimensional  distance  from  the  cone 
axis  as  shown  in  Figure  17  and  a  smooth  curve  drawn  through 
these  points.   The  value  of  the  fringe  number  in  the  region 
occupied  by  the  cone  was  taken  as  a  constant  value  equal  to 
that  at  the  cone  surface.   Prom  the  curve  so  obtained, 
values  of  the  fringe  number  at  101  equidistant  points  were 
obtained  on  each  side  of  the  cone  axis  (including  the  value 
at  the  axis)  so  that  a  total  of  201  data  values  resulted. 
This  fringe  data  was  read  into  Subroutine  RFAD  of  computer 
program  HOLOFER  and  inverted  to  obtain  the  required  density 
field. 

The  use  of  201  data  points  was  essentially  dictated  by 
the  necessity  to  be  able  to  define  the  shock  position 
accurately  and  to  have  the  fringe  distribution  in  the  region 
between  the  shock  and  the  cone  described  fairly  well.   Since 
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the  fringe  curve  from  which  these  points  were  obtained  was 
plotted  using  only  about  10  experimental  fringe  values,  the 
use  of  201  points  does  not  imply  a  higher  accuracy  in  the 
density  distribution  output  from  the  computer  program. 
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APPENDIX  R 
USE  OF  COMPUTER  PROGRAM  "HOLOFEF" 

This  computer  program  is  designed  to  invert  the  array  of 
fringe  numbers  to  obtain  the  associated  density  field  using 
the  inversion  first  proposed  by  C.  D.  Maldonado  and  described 
in  Section  III  of  this  report.   The  computer  program  can  be 
run  in  basically  3  modes  as  described  below: 

(a)  Mode  1 

This  mode  provides  the  basic  testing  capability  of 
the  program  and  uses  various  test  functions  listed  in 
Subroutine  FUNCT  in  order  to  generate  a  G  array  which  is 
then  inverted  back  to  obtain  the  original  input  function. 
Functions  other  than  those  specified  in  FUNCT  can  also  be 
read  in  on  cards  by  specifying  a  test  function  number  of 
8.0  in  the  original  list  of  42  parameters  read  in  into  the 
main  program  HOLOVERT .   In  this  case  the  main  program  first 
calls  Subroutine  FREAD  to  read  in  the  data  cards.   The  first 
card  in  this  data  deck  consists  of  two  values  indicating  the 
total  number  of  cards  to  be  read  in  and  the  'Z'  section  at 
which  inversion  is  being  performed  and  is  input  according  to 
format  89  of  the  subroutine.   Thence  follows  one  point  per 
card  according  to  format  88  representing  the  numerical  test 
function  being  input.   This  mode  was  employed  in  the  present 
case  for  inverting  the  201  points  that  were  obtained  from 
the  AGAPD  tables  by  interpolation  using  Computer  Program  1. 
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Apart  from  an  evaluation  of  the  effect  of  discontinuities 
due  to  the  shock  and  the  cone,  this  method  provided  advance 
information  on  the  magnitude  of  the  fringe  numbers  to  be 
expected  in  the  experimental  results. 

(b)  Mode  2 

This  mode  obtains  the  G  array  at  regular  intervals 
from  irregularly  spaced  values  of  the  fringe  number  function 
by  utilizing  Subroutine  SHEET.   ^he  fringe  data  may  also  be 
simulated  by  specifying  NC0DE=1  in  which  case  one  of  the 
functions  in  PUNCT  may  be  used  to  generate  the  G  array. 
This  mode  has  not  been  applied  to  the  present  investigation. 

(c)  Mode  3 

In  this  mode,  the  fringe  numbers  are  read  in  directly 
at  regularly  spaced  points  by  Subroutine  READ  which  is  called 
by  Subroutine  GARRAY.   The  various  parameters  in  the  first 
two  cards  of  Subroutine  READ  serve  to  identify  the  symmetry 
of  the  fringe  field.   The  following  parameters  have  been 
used  in  the  cases  dealt  with  in  this  report: 


Parameter      Axisymmetric  Case 

NOF        Any  value  specifying 
the  run  number. 

IMAX  201 

JMAX  1 

ISYM  2 

JSYM  101 

IMS  101 

JMS  1 

Z  The  'Z '  section  at 
which  inversion  is 
performed . 


Asymmetric  Case 

Any  value  specifying 
the  run  number 

201 

2H 

1 

1 
201 

1 

The  'Z •  section  at 
which  inversion  is 
performed. 
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Parameter 

Axisymmetric  Case 

XO 

0. 

YO 

0. 

PHISYM 

0. 

Asymmetric  Case 
0. 
0. 
0. 


Further  details  regarding  the  computer  program  are  con- 
tained in  Reference  3-   A  listing  of  the  program  is,  however, 
included  in  this  appendix  for  reference. 
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